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Abstract 

We present results of extensive computer simulations and theoretical analysis of vibrational 
phase relaxation of a nitrogen molecule along the critical isochore and also along the gas-liquid 
coexistence. The simulation includes all the different contributions [atom-atom (AA), vibration- 
rotation (VR) and resonant transfer] and their cross-correlations. Following Everitt and Skinner, 
we have included the vibrational coordinate (q) dependence of the interatomic potential. It is found 
that the latter makes an important contribution. The simulated results are in good agreement with 
the experiments. Dephasing time (r^,) and the root mean square frequency fluctuation (A) in the 
supercritical region are calculated. The principal important results arc: (a) a crossover from a 
Lorentzian-type to a Gaussian line shape is observed as the critical point is approached along 
the isochore (from above), (b) the root mean square frequency fluctuation shows nonmonotonic 
dependence on the temperature along critical isochore, (c) along the coexistence line and the critical 
isochore the temperature dependent linewidth shows a divergence-like A-shape behavior, and (d) 
the value of the critical exponents along the coexistence and along the isochore are obtained by 
fitting. It is found that the linewidths (directly proportional to the rate of vibrational phase 
relaxation) calculated from the time integral of the normal coordinate time correlation function 
[CQ(t)] are in good agreement with the known experimental results. The origin of the anomalous 
temperature dependence of linewidth can be traced to simultaneous occurrence of several factors, 
(i) the enhancement of negative cross-correlations between AA and VR contributions and (ii) the 
large density fluctuations as the critical point (CP) is approached. The former makes the decay 
faster so that local density fluctuations are probed on a femtosecond time scale. The reason for 
the negative cross-correlation between AA and VR is explored in detail. A mode coupling theory 
(MCT) analysis shows the slow decay of the enhanced density fluctuations near critical point. The 
MCT analysis demonstrates that the large enhancement of VR coupling near CP arises from the 
non-Gaussian behavior of density fluctuation and this enters through a nonzero value of the triplet 
direct correlation function. 

PACS numbers: 
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I. INTRODUCTION 



The study of vibrational phase relaxation (VPR) has been an important endeavor of 
a physical chemist /chemical physicist in the attempt to understand and quantify the in- 
teraction of a chemical bond with the surrounding solvent molecules Q, sj. As the phase 
relaxation of a bond is sensitive to the details of intermolecular interaction, it is a useful 
tool in understanding various aspects of solute-solvent interactions. The study of VPR has 
been driven by the sensitive experimental measurements of the vibrational linewidth. Thus 
the theoretical results can be verified with experiments directly. Simple models such as the 



isolated binary collisions model [,3| and hydrodynamic model Q were employed initially to 
explain the experimental results. These models led to simple expressions of temperature, 
density, and viscosity dependence of the rate. However, the agreement with experiments was 
at most tentative. It was soon realized that a major difficulty was that dephasing derives 
contributions from many sources and there may not be any unique mechanism of dephasing. 
Not only the estimation of these contributions are nontrivial, even the cross-correlations 

n 

between different pure terms are non- negligible. Thus the Kubo theory of dephasing |^ 
was straightforward to apply in this case and the calculations of the frequency modulation 
time correlation function turned out to be extremely difficult, even for simple diatoms like 
nitrogen (N2). 

Uniqueness of vibrational dephasing near the critical point is that the value of mean 
square frequency fluctuation < Ac<j^(0) > becomes large, leading to rapid decay of 
< Q{t)Q{0) >. Many factors which are responsible for this behavior are very difficult 
to understand. Near high temperature vibration-rotational (VR) coupling shows large en- 
hancement to < Auj{t)Auj{0) > including negative cross-correlation between atom-atom and 
VR coupling terms. 

Many experimental studies have been carried out on N2 usin^ vibrational Raman spec- 
troscopy as a probe. Experimental studies of Clouter et al. P, ^showed that the isotropic 
Raman line shape of simple fluid like N2 may exhibit a remarkable additional broadening 
near liquid-gas critical points {pc,Tc). They measured the Raman spectra along the triple 
point to the critical point and behavior of the line shape as the critical point is approached 
from above at a constant density. Recently Musso et al. [3] calculated the important cross- 
correlations between resonant and nonresonant dephasing mechanisms in dense liquid. They 
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observed an interesting temperature dependence A shaped linewidth (F) along the coexis- 
tence and along the critical isochore. 

In the. p,o.ee.„g stud, Oxtobv et fl fl Q showed that d.ect simulation of the 

vibrating molecules could be avoided for most cases of interest. A quantum mechanical 
perturbation theory for the vibrational motion can be used to express the dephasing rate in 
terms of auto- and cross-correlation functions of bond-force terms and its derivatives. The 
latter ones can then be calculated by molecular dynamics (MD) simulations. Oxtoby et al. 
calculated the linewidth and the motionally narrowed line shape of nitrogen near boiling 
point (77 K) 0,0]. They considered several contributions from (i) solvent-solute interaction 
force in the liquid, (ii) their derivatives, and (iii) resonant molecular vibrational interactions 



to the frequency fluctuation of mo 

Recently Gayathri et al. Q 
damental and the overtones 



ecules. 



IG [calculated the vibrational phase relaxation of the fun- 



12| | of the N-N stretch of nitrogen in pure nitrogen by 



molecular dynamics (MD) simulations. They reproduced the experimental data semiquan- 
titatively (within 40% in most cases). They have also applied the mode coupling theory 
(MOT) [iJ, to compare with the simulation results. In their calculations they have 
not include the vibrational coordinate (g) dependence of the interatomic potential and also 
ignored the cross-terms among the vibration-rotation coupling term, the atom-atom term, 
and the resonance term. More recently Everitt and Skinner studied the Raman line shape of 
nitrogen in a systematic way by including the bond length dependence of the dispersion and 
repulsive force parameters [l2|. They have also included the cross-correlation terms which 
were neglected earlier and the results for the line shift and the linewidth along the liquid- 
gas coexistence of N2 were observed to be in very good agreement with experiments. But 
calculations along the critical isochore have not been reported in their study. As mentioned 
earlier, there exists a profound experimental results in this region. 

In this work, we report results of extensive MD simulations of vibrational dephasing along 
critical isochore. We have calculated the linewidth, the line shape, and the dephasing time 
of N2 along the critical isochore and along the coexistence line. The normal coordinate time 
correlation function [CQ(t)] is calculated from the frequency fluctuation time correlation 
function for different state points along the coexistence line as well as along the isochore of 
N2. 

We have incorporated the vibrational coordinate (g) dependency of the intermolecular 
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potential and also the cross-terms. The linear expansion of the Lennard- Jones potential 
parameters on vibrational coordinate is very important to get the correct sign of line shift. 
The time integral of the diagonal and cross-terms of frequency fluctuation time correlation 
function [C(^(t)] gives the contribution to the linewidth. These cross-terms have a large effect 
on the linewidth to get a good agreement with experiment. 

The line shape calculated from the normal coordinate time correlation function shows the 
Gaussian behavior close to the critical point. Experimentally [3], it has been proved that 
the line shape remains Lorentzian for the liquid near its normal boiling point (BP). The 
increase in density fluctuations near the critical point increases the mean square frequency 
fluctuation < /\ujf >, transferring the line shape from its fast modulation, i.e., Lorentzian 
shape limit outside the critical region to a slower modulated Gaussian shape. The root mean 
square (rms) frequency fluctuation of N2 calculated along the isochore shows nonmonotonic 
behavior. However, the dephasing time (r^,) did not show any nonmonotonicity. 



II. BASIC EXPRESSIONS 



The theories of the vibrational dephasing are all based on Kubo's stochastic theory of the 
line shape. This gives a simple expression for the isotropic Raman line shape [/(tt;)] in terms 
of Fourier transform of the normal coordinate time correlation function [Cn(t)l through the 

nn 

polarizability time-correlation function as given by |4L 

/■oo 

I{uj)= I dte-xY>{iujt)[<Q{t)Q{Q)>]. (1) 
Jo 

A cumulant expansion of Eq. (Q) followed by truncation after second order gives the following 



well known expression for Cqlt) (Ref 



1)-- 



CQ{t) = < Q(t)Q(0) > 

= i?eexp(za;o^ + i < Au > t) 



X exp 



- / dt\t - t') < ALu(t')Acu(0) > 
"'0 



(2) 



where Aco'j(t) = uji{t)— < cjj > is the fluctuation of the vibrational frequency from average 
vibrational frequency. < Aco'(t)Aco'(0) > is the frequency fluctuation time correlation [Ct^(t)] 
function and ujq is the fundamental vibrational frequency of nitrogen. 
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The fluctuation in energy between the ground state and the nth quantum level of overtone 
transitions is given by 

hAu^t) = vut) - y^(o) + J2 y^iit) ■ (3) 

V^^ is the Hamiltonian matrix element of the coupling of the vibrational mode to the solvent 
bath and Ylj^iji^) represents the contribution from resonant energy transfer between two 
molecules i and j. 

The Hamiltonian for the normal mode (Q) is assumed to be of the following anharmonic 
form: 

H.^b = l^^uj^Q' + ^/g^ (4) 
2 o 

where /i is the reduced mass and / is the anharmonic force constant. The value of / is 
17.8 X 10'^ gm/cm s^. Note that Q in Eq. is not in the mass-weighted form. 

If V is the oscillator-medium interaction potential, then one finds the following expression 
for the fluctuation in overtone frequency (by using perturbation theory) Aunoit) Ref 

f^^<0it) = {Qnn - QOO) (1^) {t) + \[{Q')nn " (Q')oo] (|^) {t) 

The first two terms in the right-hand side of Eq. (0) are the atom-atom contributions and 
the third term is the resonance term. 

The vibration-rotation (VRj contribution to the broadening of the line shape is given by 
AuJnQyR = ARno/hlmre-AP Q, ^ ; whcrc AJ^it) = Pit)- < J2(0) >. However, the 
time correlation function of VR coupling term is given by 



= < A^(«)A^(0) >„„,„■« 

^)'x«J'(()J^(0)> 

"--'m ' e / 

- < J2 >2) . (6) 

J is the angular momentum and is the moment of inertia value at the equilibrium bond 
length(re). 



Auto- and cross-correlations between atom-atom forces, VR coupling and resonance terms 
have been considered in our model. The final expression for C^°*"'(t) can thus be written as 

C'f'^it) = < Auj{t)Au{0) > 

= < AuJAA{t)AuJAA{0) > 

+ < AuRs{t)AuRs{0) > 

+ < AuJAA{t)AuJvR{0) > 

+ < AuJvR{t)AujRs{0) > 

+ < AuR,{t)AujAA{0) > . (7) 

The line shape is obtained from the frequency fluctuation correlation function. F^q, 
and Fgg in Eq. (0) and the VR coupling from Eq. © are calculated separately. The main 
difference from a previous calculation |22] is that all the terms are vibrational coordinate 
dependent. This dependency comes through the bond length of the (r = r(g)) molecule. 

The Hamiltonian of homonuclear diatomic molecules can be expressed as the sum of three 
terms, 

H = H, + T(q)+U(q), (8) 

H„ is the vibrational Hamiltonian. T(q) is the total translational and rotational kinetic 
energy. U(q) is the inter-molecular potential energy, q represents the collection of vibration 
coordinates {qi}- The vibration Hamiltonian for the isolated (gas-phase) molecules is given 
by 

Here ^ represents the sum of all anharmonic oscillators for the vibrational modes of gas 
molecules. The conjugate momentum of qi is Pi for the ith molecule of the oscillator. The 
translational and rotational kinetic energy term of the oscillator can be written as 

t(q) = y + , (10) 

vh; \^2M 2I(q,) y ' ^ ^ 

where Pj and Lj are the center of mass momentum and the angular momentum for molecule 
i, respectively and 1(g) = /i(re + g)^ is the moments of inertia. We can express H = 
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T(0) + U(0) as a bath Hamiltonian and the perturbation Hamiltonian V is given by, V = 
T(q) + U(q) - T(0) - U(0). The total Hamiltonian 



H = H, +Hfe + V. 
The intermolecular potential energy can be written as 

U(q) = ^^Y1 



with 



and 



(11) 



(12) 



(13) 



(14) 



Here ijp is the unit vector along (3 atom of the vibrating molecule (zth molecule) from the 
center of mass (see Fig. (H). 

III. SIMULATION DETAILS 



We have performed a microcanonical [NVE) ensemble molecular-dynamics (MD) simu- 



lation 



23 



2J] at different state points of N2 ranging from the melting point (also the triple 



point of N2) through the boiling point and along the critical isochore (see Fig. 12)) using 
the leap-frog algorithm Q. The parameters used are giveu i,r □ Q. A system of 256 
N2 diatomic molecules were enclosed in a cubic box and periodic boundary conditions were 
used. 

The thermodynamic state of the system is expressed in terms of the reduced number 
density of p* = pa^ (Ref. 2^) and a reduced temperature of T* = ksT/e. a is the diameter 
of the molecule and e is the interaction parameter (seeH}. The unit of temperature is e/fc^ 
(K), where fc^ is the Boltzmann constant. Cheung and Powels j27| had earlier studied liquid 
N2 at different state points using MD simulations. Most of the thermodynamic state points 
chosen for the work presented here have been taken from their study. We have done few 
simulations with a system of 512 nitrogen molecules to check the system size dependency. 

Fig. 121 gives a schematic view of the phase diagram 3|. The arrowed line points out 
that along the critical isochore, is approached from above. For nitrogen the triple point 
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non-vibrating molecule 
P th atom 



vibrating molecule 



CLth atom 




^th atom 



ongin 



FIG. 1: The schematic diagram illustrating molecular interactions between the nonvibrating (zth 
molecule) and vibrating (jth molecule) homonuclear diatomic (N2). 



corresponds to that given by {T^.pD — (1.7212,0.6964) and the critical point, {T*,pl) — 
(3.3592,0.35). 

For intermolecular potential-energy (Vjj) between two molecules i and j, the following 
site-site Lennard- Jones type is employed as given below, 

1,2 

V^=EV(W)- (15) 
Here V(rjQ,j^) is the Lennard-Jones atom-atom potential defined as 



(16) 



Vibrational coordinate dependence of e and a has been incorporated following Everitt and 
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TABLE I: Parameters for N2. 




' 01 


1.094 




tlkK 

1 


37.3 




a/°A 


3.31 


Spectroscopic Constants 


M/amu 


28.0 






2358.57 


Polarizability parameters 


7/1-^ 


0.62 




5/1-1 


-0.063 



[3.3592] 



[1.7212] 




[ 0.6964 ] 



(reduced density) 



FIG. 2: Phase diagram of a typical substance, showing boundaries between solid (S), liquid (L) 
and vapor (G) or fluid (F) phases. This is only the projection of the p* — T* plane. The parameters 
(T*, T/, p*, and Pf) are given for nitrogen. We did the simulation along the critical isochore of 
nitrogen which is indicated by an arrow. 
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Skinner 



m 



e{q) = e[l + 27g]; 

a{q) = a[l + 25q]. (17) 



For homonuclear diatomic-like nitrogen, 



= (18) 

^iajP — 2 

= (19) 

and 

~ e(l + + ^qj + 2-f^qiq^), (20) 

aij =a{l + 6qi + 6qj). (21) 
Now the LJ potential takes the form as below, 

2 

Q,/3 = l 

x{i^Y\l + 5q, + 5q,Y' 

-i^m + Sq. + 6q,n]. (22) 

IV. SIMULATION RESULTS, COMPARISON WITH EXPERIMENT AND DIS- 
CUSSION 

A. Along the coexistence line 

The frequency-modulation time correlation function [C^j (/:)], the dephasing linewidth, 
and the line shift j^l are all obtained for several thermodynamic state points ranging from 
triple point to critical point. 

We have calculated the line shift Au = {{uji) — ujq)/27cc as a function of density for 
nitrogen. The magnitude of the line shift increases with density as clearly seen in Fig. Efa) 
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TABLE II: Simulated values of Dephasing time at different temperatures. 



Temperature(K) 




r„ fps) 




Simulation 




Expt 


125.3 


12.40 




14.5 


140.0 


19.06 




23.3 


149.3 


25.13 




27.8 


170.0 


26.62 




28.6 


186.5 


26.75 




26.7 



and the results are in good agreement with the experiment. The vibrational coordinate 
dependence is important to get the correct sign of the line shift. 

Along the coexistence line the root mean square frequency fluctuation (A) increases as 
the critical point is approached. Fig. El^b) shows that A decreases linearly with density with 
a slope of -0.3624. 

B. Along the critical isochore 

We have plotted C^^ii) against time for three different temperatures along the critical 
isochore in Fig. 0] The lines with open circles show Ci^if) which decays faster than C^^ify 
represented by the simple line. 

The temperature dependence of the dephasing time (xy) calculated by integration of C^(t) 
[Eq. Cj is comparable with experimental results (see ITT|. 

In Fig. El a simple line represents the normal coordinate time correlation function without 
including q dependency of L J parameters whereas the dashed line represents the same results 
with q dependent L J parameters. The former one has a larger correlation time than the latter 
one. The Cqii) clearly shows the importance of q dependence. 

Fig. ini displays the root mean square frequency fluctuation A [= AcJ^(O)] calculated 
from unnormalized C^(0) at different temperatures along the isochore. The calculated A 
shows a nonmonotonic behavior with temperature. A initially increases with temperature 
and shows a maximum near critical temperature indicated by an arrow. A sharp decrease is 
observed for the temperature greater than T^. The increase in density fluctuation near the 
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0.3 0.4 0.5 0.6 0.7 0.8 0.9 
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* 

p 



FIG. 3: The Raman line shift {Au) is plotted against density (p*) (a) along the coexistence line 
for N2. The critical density is indicated by the arrow. The rms frequency fluctuation A is plotted 



against reduced density (5) along the coexistence line for N2. Simulation results have been fitted 
with the formula y = ax + h with a = —0.3624 and h = 0.3064. 
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0.4 

time (ps) 

FIG. 4: The frequency fluctuation time correlation functions are plotted against time 

(t) at different temperatures (a) T = 45.0 K, (b) T = 140.0 K, and (c) T = 300.0 K along the 
critical isochore. In all the figures above, circles with a simple line show Ci^{t) for q dependent 
interaction potential, whereas the simple solid line shows the q independent interaction potential. 

critical point is responsible for this nonmonotonic behavior of rms frequency fluctuations. 

Linewidth shows an interesting A-shaped feature when it is plotted against temperature 
(see IIII|) for the two different branches of calculation for N2 as shown in Fig. [7| This figure 
is very similar to the one observed in experiment (see Fig. |3]of Ref 0). It is interesting 
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-0.2 

20 40 60 80 100 120 140 150 

time (ps) 

FIG. 5: The normal coordinate time correlation functions [C^Q^t)] are plotted against time {t) 
at three temperatures along the critical isochore (a) T = 110.8 K, (b) T = 125.3 K, and (c) T = 
140.0 K. 

to note the sharp rise in the dephasing rate as the critical point is approached. There are 
six main contributions, (a) density [C'^(/:)], (b) VR coupling [C^^{t) ], (c) resonance 
[C^'{t)], (d) density-VR coupling [C^"^^(t)], (e) density-resonance [C^^'^^t)], and 
(f) VR-resonance [C^~^^{t)] which are responsible for the sharp rise in total linewidth near 
the critical point. These are the time integrals of diagonal and cross-terms in the frequency 
fluctuation time correlation function. We shall come back to this point a bit later. 

A crossover from Lorentzian-type line shape to Gaussian line shape is found to take place 
when there is a large separation in the time scales of decay of C^(t) and Cglt) 0, 1^ cease to 
exist and the two time correlation functions begin to overlap. We have calculated Ci^{t) and 
Cgit) at three temperatures near the critical point (see Fig. |HI). The decay of C^(t) becomes 
significantly faster, reducing the gap of decay between the two correlation functions. Indeed, 
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0.2 




T(K) 

FIG. 6: The rms frequency fluctuation A vs temperature (T) along the critical isochore for N2. 
We fitted all simulation points with Gaussian. The fitting formula is y = a[e~^^^~'^^'^], where 
a = 0.184904,6 = 0.00010813, andc = 118.874. 

the computed line shape becomes Gaussian near the critical point but otherwise remains 
Lorentzian-type both above and below the critical temperature. Note that the frequency 
modulation time correlation function decays fully in about 200 fs. 

To understand the origin of this critical behavior, we have carefully analyzed each one of 
the six terms which consist of three autocorrelations and three cross-terms between density, 
vibration-rotation coupling, and resonance terms which have been mentioned earlier. 

Fig.Elshows the time dependence of the four dominating terms , the two autocorrelations 
[Fig. IHK^a) and Fig. El(b)] and two cross-correlations (Fig. Efc) and Fig. Efd)). The 
decompositions of line shift, linewidth, and the temperature dependent quantities are fully 
dependent on the contributions which come from all six terms. As we approach the critical 
point along the critical isochore, the magnitude of contributions from different terms are 
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TABLE III: Temperature dependence A shaped Unewidth along the coexistence and critical iso- 



chore. 


Temperature(K) 


Linewidth(Ghz) 


64.2 


2.764 


71.6 


2.802 


75.7 


2.809 


89.5 


2.812 


110.8 


2.812 


120.0 


6.965 


125.3 


16.12 


135.0 


12.38 


140.0 


10.49 


149.2 


7.96 


170.0 


7.51 


186.5 


7.47 



increased. Results for only four different temperatures including the critical point are shown 
in Fig. ini All the remaining terms, resonance-resonance and density-resonance, are found to 
be unimportant in comparison with these four terms. 

The time dependence relative contribution Xij{t) = j^dt'{t — t')C^J{t') of the density 
[Fig. 1 10b ] and vibration-rotation coupling term (Fig. ITUb) are plotted for seven state 
points along the critical isochore, which are found to be dominant near the critical point. 
The sharp rise in the value of the integrand as the critical temperature is approached and 
the fall when it is crossed. We have found that both these contributions at the critical point 
are distinct compared to the other state points. Thus the rise and fall of the dephasing rate 
arise partly from the rise and fall in the density and the vibration-rotation terms. We have 
calculated the slopes of the relative contributions by the linear fitting in the long time at a 
particular temperature T = 125.3 K. The values of the slope are 5.4 x 10~^, 9.5 x 10~^ and 
—6.8 X 10~^ for density-density, VR-VR coupling, and density- VR terms, respectively. The 
major contribution comes from the VR coupling term among the four dominating terms. 
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25 




60 80 100 120 140 160 180 200 



T(K) 

FIG. 7: The lambda (A) shaped hnewidth (F) of nitrogen along the coexistence line (closed square) 
and the critical isochore (closed circle). The experimental results for linewidth along the coexistence 
line and the critical isochore (open circle) reported by Musso et al. in Ref. P| are also shown. The 
critical point is indicated by an arrow on the abscissa. 

V. DYNAMICAL HETEROGENEITIES NEAR THE CRITICAL POINT. 



We have investigated the presence of dynamical heterogeneities in the fluid to further 

This can be 



explore the origin of these anomalous critical temperature effects Isol 
quantifled by the well-known non-Gaussian parameter a{t). It is deflned as 

'3\ < r\t) > ^ 



a{t) 



(23) 



5 J <r2(t) >2 

where < Ar(t)2 > is the mean squared displacement and < Ar(t)'^ > the mean quartic 
displacement of the center of mass of the nitrogen molecule. It can only approach zero (and 
hence Gaussian behavior) at a time scale larger than the time scale required for individual 
particles to sample their complete kinetic environments. The function a{t) is large near the 
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1?^ 




time (ps) 

1* ' ' 




time (ps) 

la , , 




time (ps) 

FIG. 8: The frequency fluctuation time correlation function, Cuj{t), and the normal coordinate time 
correlation function, CQ{t), are plotted at (a) 186.5 K, and (b) 140.0 K (along critical isochore), 
and (c) the Cuj{t) and Cqlt) are plotted along the coexistence line at temperature 75.7 K. 

critical point at times 0.5-5 ps as observed in Fig. ^2 which indicates the presence of long 
lived heterogeneities near T^. 
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FIG. 9: The decopositions of unnormalized frequency fluctuation time correlation functions for 
N2 into (a) density-density, (b) VR-VR coupling, (c) density- VR and (d) VR-resonance at four 
temperatures along the critical isochore. The density-resonance and resonance-resonance terms are 
relatively small and are not shown here. 
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0.8 1 1.2 

time (ps) 




0.8 1 1.2 

time (ps) 



FIG. 10: The time dependence of the relative contributions of the (a) density-density, X^gnit), and 
(b) VR coupling, XvR{t), terms at different temperatures along the critical isochore. 
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FIG. 11: The non Gaussian parmeter a{t) is plotted against time (t) at three different temperatures 
along the critical isochore. 

VI. NEGATIVE CROSS-CORRELATION BETWEEN THE DENSITY AND VR 
COUPLING TERM 

The role of cross-terms (which can be negative) are extremely important for the decay 
of C^{t) which occurs in the femtosecond time scale along the critical isochore. We have 
calculated the cross-terms among density, VR coupling, and resonance term. The cross- 
terms of VR coupling with density and resonance are found to be negative [see Figs. ISl(c) 
and|ni(d)] . As mentioned earlier the major contributions to the frequency fluctuation come 
from the density as well as from the VR coupling term. One of the main reasons for ultrafast 
decay of C^j is the cancellations of negative cross-terms from the total contribution. 

Now the question is why are these cross-terms are negative? The distribution of density 
dependent and the VR coupling dependent part of the total frequency fluctuation along the 
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FIG. 12: The distribution of fluctuation in (a) the density term, and (b) the VR coupHng contri- 
bution to total fluctuation in the frequency at the temperature 110.8 K. 
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isochore are shown in Figs. IT^ a) and lT^ b). respectively. The peak of those distributions is 
near zero 5ujvr which proves that the homogeneity condition is satisfied along the isochore. 
It is evident that the distribution of 5ujvb. is always positive while the distribution of 5up 
showed the long negative tail. The distribution of the product of density and VR coupling 
terms are plotted for different temperatures in Fig.^J It is clearly seen that this distribution 
is negative. A detailed analysis of the origin of the negative sign of the cross-correlation 
between AA and VR coupling terms is given in Appendix B. On the other hand Aujp is 
directly proportional to the force acting on the bond of the diatoms. If force F is large, the 
velocity will be large and it reflects that the will be decreased when F is increased and 
vice versa. This is the origin of anticorrelation between the density and VR coupling terms 
in C^{t). 

The distribution of correlation between the density-density term and the VR-VR cou- 
pling term in the frequency fluctuation time correlation function are shown in Figs. IT^ a) 
and ITlT b). respectively. Both distributions are positive and the homogeneity condition of 
frequency fluctuations is satisfled in both cases. 



VII. MODE COUPLING THEORY ANALYSIS. 



We can use mode coupling theory to obtain an expression for the atom-atom contribution 
to the frequency modulation time correlation function. The main steps have already been 
discussed by Gayathri et al. and need not be repeated here. In brief, MCT gives the fol- 
owing expression for the density dependent frequency modulation time correlation function 



33[ considering that the number density is the only relevant slow variable in dephasing. 



(Ao;,(0) Aa;,(t)) = H edkF,{k, t)c\k)F{k, t). (24) 

Where c{k) is the Fourier transform of the two particle direct correlation function. The main 
contribution is derived from the long wavelength (that is small k) region near the critical 
point (CP). In the small k limit the hydrodynamic expressions for two-point correlation 
functions are given by 

Fs{k,t) = e'""^^'', (25a) 
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FIG. 13: The distribution of correlation between density- VR coupling at three temperatures, (a) 
110.8 K, (b) 125.3 K and (c) 140.0 K, respectively. Note that this figure is entirely different from 
Fig.d 



F{k,t) = S{k)e 



(25b) 



where Fs{k,t) is the self-intermediate scattering function and F{k,t) is the intermediate 
scattering function. Here Dg is the self-diffusion coefficient, S{k) is the static structure 
factor, and Dt is the thermal diffusivity. Thus {Aujp{0)Aujp{t)) ~ S{k o)e-(^-+^^)^'*. 
Near the CP, S{k — > 0) becomes very large ( as compressibility diverges at T = Tc ), 
leading towards a Gaussian behavior for line shape. Dt also undergoes a slowdown near Tc. 
However, a limitation of the above analysis is the absence of the VR term which contributes 
significantly and may mask some of the critical effects. 
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three temperatures 110.8 K, 125.3 K and 140.0 K respectively. 



All these lead to a large value of < Au;^(0) >, which may lead to a Levy distribution 
from time dependence of {Q{0)Q{t)) as discussed earlier by Mukamel et al. 

H- At high 

temperature, the latter dominates over the density term. We have calculated the line shape 
and fitted with the Fourier transform of Levy distribution ( see Fig. [T3j) . Clearly as is 
approached, the low k fluctuations become more important and, if the high k contributions 
are sufficiently small, they will dominate the line shape. We would like to mention that 
the line shape becomes Gaussian (see the Fig. ITH|l instead of the Fourier transform of Levy 
distribution near the critical point. The small k picture described here is quite from the usual 
collisional broadening picture of dephasing. A complete Lorentzian behavior is predicted 
only in the low temperature liquid phase. Interestingly, the predicted divergence of < 
Aco'^(O) > very close to Tc enhances the rate of dephasing and this shifts the decay of 
< Q{t)Q{Qi) > to short times, which gives rise to the Gaussian behavior. 
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FIG. 15: The behavior of hne shape at the temperature 123.5 K. 

We next use MCT to analyze the cross-correlation between VR coupling and the atom- 
atom (AcUp) term. According to density functional theory (DFT) the effective potential 

energy (Vg//) of the mean field force [F = — VVe//(r,t)] can be written in terms of the 
two-particle correlation function c(r) without including orientation(f2) as 

/5K//(r, t) = -J dr'c{T - r')6p{r', t). (26) 

Assuming that the density term is approximated by the isotropic limit, we can therefore 

write the atom-atom terms as 

A6Up(r, t) = -keT/h J dr'c{r - r', t)Sp{r', t). (27) 
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The calculation of the VR term is a bit more complicated. The angular momentum J{t) at 
time t can be expressed in terms of torque N(t') at earlier times {f > t), 



Jit) = J(0) + J dt'N{t'). (28a) 

From Eq. EHlof the Appendix B, the VR coupling can be written as 

AuvniT, n, t) = A3{v, ft, t) ■ J(r, Q, t) = A[J^{r, n,0) + 2 j dt' J{r, U, 0) • N(r, n, t') 

+ j dt' j dt"N(r, O, t') ■ N(r, O, t")](28b) 

Here A is a constant (see Appendix B) and (5p(r', f2', t') is the fluctuation in the number 
density. 

The expression for the cross-correlation can be expressed as 

(Acuyfi(r, O, t) Aa;,(r, 0)) = -A{kBT/h) [ j dv'c{v - r') {j\v, O, 0)(5p(r', 0)> 

+2 j dt' j dr'c(r-r') (J(r, O, 0) ■ N(^^, 0)) 

+ j dt' j dt" j dr'c{r - r') (N(rj, t') ■ N(0, t")6p{r', 0))]. (29) 

By using DFT again, the torque (N) can expressed in terms of density functions as 

N(r, n, t) = Vnj dr'dn'c{r - r', Q, 0')5p(r', Q', t). (30) 

Combining the above equations, the final expression cross-correlation can be written as, 

< AuJvRir, n, t)Au,ir, 0) >= -AiksT/h) [ J dv'cij - r') {j\v, O, 0) V(r', 0)> 

+2 j dt' j dr' J dr"dn'c{r - r') Vn' ■ c(r - r", fJ, ft') (J(r, fi, 0)5p{r", n', t')5p{r', 0)) 

+ j dt' J dt" J dr"dn' j dv'"dQ," j dv'c{r - r') Vn' ■ Vn"c(r - r", O, n')c{r - r'", O, Vt") 

x(5p(r",O',t')5p(r'",O",t")V(r',0))] 

= I + II + II|31) 

We now analyze Eq. (|HT|) term by term. Since both J and 5p are uncoupled at the same 
time, the contribution of I is zero. The second term II in Eq. fj31|) consists of two density 
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terms with an angular momentum term. Using cumulant expansion we can show that the 
contribution of II is also equal to zero. 

To evaluate the third term (III), we factorize the translational and angular variables in 
the two-particle direct correlation function (DCF) and we can write 

c(r - r', ft') = ce{ft, ^2')co(|r - r'|). (32) 

Where Co(|r — r'|) is the isotropic part of the direct correlation function. When the temper- 
ature T — i> Tc, Co Po ^, which is finite. The angular part of the direct correlation function 
co{fl, fl') can be expanded in terms of the spherical harmonics. 

ce{n,n')= J2 ('hhmYi,,m{n)Yi,,^{n'). (33) 

Zl ,l2 ,171 

If we assume that the density effects are included in the spatial part co(|r — r'|) then the final 
expression of the angular part of the direct correlation function will be finite and nonzero. 
If one assumes that the 6p are Gaussian random variables, then the value of the triplet 
correlation function in III is zero, by cumulant expansion. However, the distribution of 
the fluctuation in number density clearly shows a non-Gaussian behavior, arising from the 
long lived heterogeneity near critical point (see Fig. 11). We, therefore, expect a nonzero 
value of the cross-correlation between VR coupling and atom-atom terms from the triplet 
correlation function in Eq. ()31|). Moreover, the amplitude of the cross-correlation is predicted 
to be small at all temperatures away from Tc but it is predicted to become significant as Tc 
is approached, where the density fluctuation becomes non-Gaussian. 

We have calculated the term AR in Eq. EH (see Appendix B) and found it to be always 
positive. It is obvious that > 0. So we can conclude that 6ujvr > [see Fig. 126]. 
Figure 12(a) shows that most of the distribution of fluctuation in 6ujp are in the negative 
direction. On the other hand the distribution of fluctuation in Suvr is in the positive 
direction [see Fig. 126]. Here, Aup = {nh{-f)/2ij,^uj^) F^q + {nh/2fiuJo) F^q. So the 
distribution of cross-correlation between density and VR coupling terms must be negative 
(see Fig. 13). 

The above analysis MOT demonstrates that the large enhancement of vibration-rotation 
coupling near the gas-liquid critical point arises from the non-Gaussian behavior of density 
fluctuation and this enters through a nonzero value of the triplet direct correlation function 
[Eq. (EH)]. 
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VIII. DIVERGENCE OF RAMAN LINEWIDTH NEAR THE CP 



Along the coexistence line and the critical isochore, we have found that the temperature 
dependence of the linewidth (F) is singular near the critical point. The divergence-like 
rise near is fitted to a form (T — T^)"^'. From the fitting with both experimental and 
theoretical data we have found /?' = 0.386 along the coexistence line and /3' = 0.207 along 
the critical isochore. 

Above Tc, the fluctuations are small which cause the resonance lines to be homoge- 
neously broadened and the line shape is Lorentzian. On approaching T^, the fluctuations 
35I l3q become large and their correlation time also increases, as a consequence the line 
shape becomes Gaussian. 

The order parameter for the liquid-gas critical point is (pl — pa)- This goes to zero as 



(Tc - T)-^ where ~ |. Mukame: 
lines by using mode-coupling theory 



et al. (3J] have studied the broadening of spectral 
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38| near a liquid-gas critical point. They assumed 
Lorentzian form of the line and found that the linewidth increases with e~^, where e = 
|(T — Tc)/Tc|. They have also found that the value of critical exponent (s) changes from 
0.607 for large e to for small e. That is, the singularity becomes weaker as is approached. 
It is likely that simulations are not sufficiently close to Tc because simulations are finite sized. 

The reason for the relative success of a small system in reproducing experimental behavior 
is that dephasing probes only local static fiuctuations. Thus long range static and dynamic 
correlations important in specific heat, compressibility, or light scattering are not relevant in 
vibrational dephasing. Second, while our simulated system can certainly capture the increase 
in fiuctuation, it can never capture the long range fiuctuations. But it captures enough to 
reproduce many of the features. Note that even a small system is capable of exhibiting 
large fiuctuations near the gas-liquid coexistence/critical point. Put bluntly, truly large 
scale critical fiuctuations are neither required nor refiected for dephasing simply because 
dephasing is a local process. 

IX. CONCLUSION 

In this article extensive MD simulations of vibrational phase relaxation of nitrogen have 
been presented. The simulations reported here seem to reproduce many of the anomalies 
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observed in experiments. It shows (a) the origin of enhanced negative cross-correlations, (b) a 
crossover from a Lorentzian-type to a Gaussian hne shape as the critical point is approached, 
(c) the non-monotonic dependence of rms frequency fluctuation on temperature along the 
critical isochore, (d) a lambda shaped dependence of linewidth on temperature, and (e) a 
near divergence of linewidth near the critical point. 
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XI. APPENDIX A: DERIVATIVE OF LENNARD-JONES POTENTIAL 



The derivative of the LJ potential can be calculated as 



dqi 



4e J2 [{1 + 7(1 + Sij) + 2j\qj + q^j)} x { ^ (1 + % + Sq^)'' 
(1 + H + + {1 + l{<li + (Ij) + 27'?i?j} 



X { 



il + Sqi + SqjY 



dr-ia^jfi XTiaj^J dq, 



(— )°X6(1 



+ 5 + 55ij){l + 5 + 5qi + 5qjf}l 



(34) 



[where 



12 



12 



13 



and 



d 



dr. 



ictjl3 \ 'iaj0 
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4e ^ [{1 + 7(1 + bij) + 272(5, + q^bij)} X { ( ^ 1 (1 + bq^ + 

- (l + 5?i + %f} + {l + 7(?i + ?,) + 2725,5,} 

X {7^f^y'%^(i+%+<^..)" 

+ (^) (l + (5 + 5(5i,) X 12(l + % + (55,)^^ 

- ( ^ ) X 6(l + 5 + M,,) X (l + 5 + 55, + 55,)^}]. (35) 

XII. APPENDIX B: FREQUENCY MODULATION FROM VIBRATION- 
ROTATION COUPLING 



The vibration- rotation centrifugal coupling term is Eyib^rot = where /(= 2/ir2) is 
the moment of inertia and fi is the reduced mass of a diatomic molecule. This VR term is 
important only when I depends on vibrational coordinate (g) through the bond length, i.e., 
r{q) — rQ + Ar{q). Expanding this as a Taylor series about the equilibrium bond length ro, 
the Eyib^rot can be rewritten as 

J2 



vib—rot — 



2 



2//(ro + Ar) 

1 + — 



J2 



^_2Ar^3/Ar- 



(36) 



27o 

The contribution of VR coupling to the total Ci_^{t) is given by C^^ as 

Cj^(i) = {Au;vR(t)Au;vRm (37) 

where 

Aa;vi^(t) = [Eu{t)-Eoo{t)]/h 
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where AR = (Qn - Qoo) - n\Ar\n > and = {n\{Arf\n) are the 

expectation values of bond length displacement and its square and A J^(t) = J^(t) — ( J^(0)). 
We used (J^) = 2/oA;bT. 
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